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^ , Abstract 

Ph' We consider the compatibility of DAMA/LIBRA, CoGeNT, XENONIO and 

Q . XENONIOO results for spin-independent (SI) dark matter Weakly Interacting Mas- 

+i , sive Particles (WIMPs), particularly at low masses (~10 GeV), assuming a standard 

I dark matter halo. The XENON bounds depend on the scintillation efficiency factor 

£cff for which there is considerable uncertainty. Thus we consider various extrapo- 
CN ■ lations for Ces at low energy. With the Ces measurements we consider, XENONIOO 

^ ■ results are found to be insensitive to the low energy extrapolation. We find the 

, strongest bounds are from XENONIO, rather than XENONIOO, due to the lower 

energy threshold. For reasonable choices of CeS and for the case of SI elastic scat- 
tering, XENONIO is incompatible with the DAMA/LIBRA 3cr region and severely 
• constrains the 7-12 GeV WIMP mass region of interest published by the CoGeNT 

■ collaboration. 



H ■ I. INTRODUCTION 

The nature of the dark matter that comprises a quarter of the Universe is one of the 
big unanswered questions in astrophysics and particle physics. Perhaps the best motivated 
candidates are Weakly Interacting Massive Particles (WIMPs) which have weakly interacting 
cross sections and masses in the GeV-10 TeV range. In recent months, there have been new 
data releases from many experiments that have engendered a great deal of excitement. Of 
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particular interest is a low mass region ~ 10 GeV which at first sight seems to be compatible 
with a number of different experiments. The goal of this paper is to examine some of the 
issues regarding the question of such compatibility. In this paper, we will restrict ourselves 
here to considering WIMPs with spin-independent (SI) interactions. Spin-dependent and 
mixed couplings will be examined in a future work. 

It was initially the DAMA/NAI experiment fv\, looking for annual modulation j2|, [sj of a 
WIMP signal, that found such a possible low mass region. WIMPs with SI interactions in 
the mass range 5-9 GeV were found to be compatible with the DAMA/Nal results and all 
negative results from other searches that existed at the time 0, [sj . The situation changed 
after the publication of the first DAMA/LIBRA results {oj (see e.g. Ref. 0] and reference 
therein). For SI interactions, Ref. |7j] found that the best fit DAMA regions were ruled out 
to the 3a C.L. But Ref. also found that for WIMP masses of ~8 GeV, some parameters 
outside these regions still yielded a moderately reasonable fit to the DAMA data and were 
compatible with all 90% C.L. upper limits from negative searches, when ion channeling in 
the DAMA experiment as understood at the time was included (see Section [II Bl below). 
The strongest bounds at the time came from CDMS jsf and XENONIO |9i]. Since then 
many new data sets have been released and a reexamination of the light WIMP region 
is now necessary. We will focus on the following three factors: a possible dark matter 
signal for 7-12 GeV WIMP's found by the CoGeNT collaboration Il0| : the existence of new 



better upper limits, in particular the bounds set by the XENONIO [Ul and XENONIOO [12 
collaborations; and the recognition that the effect of channeling in Nal(Tl) crystal is less 
important than previously assumed [l^ (see Section Hi Bp . 

CDMS [1^ has released its full data set, with tighter bounds and two unexplained events 
at low energy that may be compatible with background. CDMS constraints will be included 
but are not a focus of this paper. Our focus is on elastic scattering of spin-independent 
WIMPs from a standard Maxwellian halo; recent examinations comparing experimental 



studies in this case include 15-22 



Just prior to our paper being submitted, a revised version of Ref. |2J| appeared that 
added an examination of XENONIOO and Ces in the context of scalar WIMPs. The halo 
model parameters and £eff models used in that paper differ somewhat from ours, but the 
results are quahtatively similar. 

The interpretation of the XENONIO and XENONIOO results requires the ability to re- 
liably reconstruct the nuclear recoil energy from the observed signal. This depends on the 
scintillation efficiency factor Ces for which there is considerable uncertainty at low energies 
(see Section Hi A|) . A large part of this paper is devoted to examining the Ces dependence 
of the two XENON constraints. 

In this paper we focus on comparing the following experimental results: the combined 
modulation signal fl^ as well as the total rate ^ from DAMA/Nal and DAMA/LIBRA; 
the combined CDMS 5-tower results 0, [l3l; the recent first resuhs from XENONIOO [Ti]; 
and the older but lower threshold XENONIO reanalysis results [ll|. Constraints for CDMS, 
XENONIO, and XENONIOO are determined using the Maximum Gap method^ 2^, while 



^ For zero observed events (as in the case of XENONIOO), the Maximum Gap (MG) method provides an 
identical constraint as that produced by a Poisson hmit based on the total number of events. When events 
are observed, the MG method provides better constraints than the Poisson case as the former takes the 
energy spectrum into account, whereas the latter does not. The use of MG is of particular importance for 
XENONIO as the 13 observed events have energies that are inconsistent with the spectrum expected for 
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the parameters compatible with DAMA are determined via the goodness-of-fit of their ob- 
served modulation signal with the theoretically expected signal. Details of these statistical 
analyses may be found in Ref. 

For the two XENON experiments, we assume the energy resolution is primarily limited 
by a Poisson distribution in the small number of photoelectrons (PE) expected at low re- 
coil energies. Interactions in the liquid Xenon comprising the XENONIO and XENONIOO 
detectors give rise to a prompt scintillation signal (51) followed by a delayed secondary 
scintillation signal (5*2). The quantities 5*1 and S2 are discussed in the following section 
and the various thresholds and data cuts are described in e.g. Ref. The efficiencies for 
XENONIO and XENONIOO are taken from Refs. 11| and 12[, respectively. However, for 



XENONIO, the SI peak finding efficiency factor rjsi, which was not included in Ref. [11 
(where it was not particularly relevant), must also be taken into account. For this rjsi factor, 
we take the more conservative of the two cases found in Ref. [25j . 

A second issue must also be taken into account. As explained in detail below, as the 
recoil energy decreases, so do the average 5*1 and 5*2 signals. At low enough energies, a 
sizable fraction of the events may fail to produce enough 5*2 signal and/or to fall in the 
proper log(S'2/S'l) range (the nuclear recoil band cut), even if a fluctuation gives a high 
enough 5*1. At higher recoil energies, the relative size of the fluctuations get smaller and the 
average 5*2 is too high for any signiflcant fraction of the 5*2 fluctuations to fall below the 5*2 
threshold. However, somewhere below an 5*1 average of 1 photoelectron (PE) in XENONIO, 
an event would need not only an upward fluctuation in 5*1, but an upward fluctuation in 5*2 
to pass both thresholds and it would require particular ranges of fluctuations to fall within 
the required range for log(S'2/S'l). Thus, to avoid issues with the 5*2 threshold and nuclear 
recoil band cuts, we ignore recoil energies that give on average less than 1 PE in the 5*1 
signal. We flnd that, with Poisson fluctuations included, low energy recoils in XENONIO 
tend to pass the various thresholds and cuts at a higher rate than predicted by the efficiencies 
we use; these efficiencies can thus be considered conservative over the range at which they 
are applied. Additional events passing these various cuts can arise from low energy recoils 
that give an average SI signal (5*1) that falls below our imposed cutoff of {SI) > 1.0 PE. 
Accounting for these events will strengthen the XENONIO constraints at low WIMP masses. 
Our (5*1) cutoff is simply due to the efficiency being unknown for these low energy recoils 
at this time. We are examining these low energy recoil efficiencies for a future work, though 
a comprehensive treatment of these efficiencies in XENONIO and XENONIOO has appeared 



in the meantime 26 



In addition to constraints for the above experiments, we also show the 7-12 GeV WIMP 
mass region suggested by CoGeNT as an explanation for excess events seen at low energies in 
their detector [lO[. We perform our own statistical analyses of all data sets except CoGeNT, 
for which we simply use their published region. 

To allow for direct comparison, all the other experimental constraints are determined us- 
ing the same 600 km/s galactic halo escape velocity as used by CoGeNT in their analysis (27j . 
This escape velocity falls within the 90% confidence interval of 498 km/s to 608 km/s that 
was determined by a recent analysis of high velocity stars [28| , but is somewhat above the 
median likelihood of 544 km/s that was found in that analysis. Similarly, we use throughout 
the value 220 km/sec (used by CoGeNT) for the rotation velocity of the Galactic disk in the 
vicinity of the Sun. We note that new measurements suggest that the rotation velocity might 



a light elastically scattering WIMP. 
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be higher, 254 ± 16 km/s [29|, the effect of which is to shift the best fit in all experiments 
to slightly lower WIMP masses [sO*] . 

Comparison with Version One of this paper: We have made two substantial 
changes, as discussed in the previous paragraphs. First, we have included the effects of 
risi, the 51 peak finding efficiency factor. With this factor included, the XENONIO con- 
straints shift upwards in scattering cross-section by at most ~ xl.7 at any given WIMP 
mass, not by several orders of magnitude (as was speculated in Ref. [SlJ)- Second, in the 
current version we do not include recoil energies which yield on average less than 1.0 PE in 
the 5*1 signal. It is this second issue that is by far the bigger effect. A more complete dis- 
cussion can be found in the Appendix. In addition, in the Appendix we respond to various 
critiques of our work. 



II. EXPERIMENTAL ISSUES: ^eff AND CHANNELING 

In this section we discuss two important experimental issues. First, the scintillation factor 
in XENON C^s is extremely important in interpreting results and yet is not well known. 
Second, the channeling effect in DAMA/LIBRA, again not well known, may change the 
location of the regions in WIMP parameter space that are compatible with the data. 



A. The £eff Scintillation Efficiency Factor in XENON 

The interpretation of the XENONIO and XENONIOO results requires the ability to re- 
liably reconstruct the nuclear recoil energy from the observed signal. This reconstruction 
depends on the scintillation efficiency factor C^s for which there is considerable uncertainty 
at low energies. Here we discuss this factor and present three models for C^s at low energy. 

Interactions in the liquid Xenon comprising the XENONIO and XENONIOO detectors 
give rise to a prompt scintillation signal, 5*1, followed by a delayed secondary scintillation 
signal, S2. The SI signal arises from a rapid relaxation of excited Xenon states produced 
as a result of the interaction. The 5*2 signal arises from ionized electrons also produced 
in the interaction; these drift through the liquid xenon under an applied electric field, but 
once they reach the liquid surface they are extracted into a xenon gas phase where they 
emit proportional scintillation light. The drift time of the electrons causes this secondary 
scintillation (the 5*2 signal) to be observed later than the SI signal, allowing both scintillation 
signals to be measured separately. The 5*1 signal can be used to determine the energy of 
the interaction, while the combination of both signals allows discrimination between nuclear 
recoil events (possibly WIMP interactions) and electron recoil events (necessarily background 
interactions). The ratio of S2 to SI is much higher in the case of electron recoils than in 
the case of nuclear recoils. 

Interpretation of the XENON results requires the ability to reliably reconstruct the nu- 
clear recoil energy Enr from the observed 5*1 signal. Calibration of the nuclear recoil energy 
dependence of 5*1 often involves gauging the detector's response to electron recoils at higher 
energies; parts of the detector's response {e.g. the fraction of scintillation photons that yield 
photoelectrons (PE) in the photodetectors) are more easily determined in this case than 
with nuclear recoils at lower energies. Taking 5"! to be normalized to the number of PE, 5*1 
and Enr are related by an equation involving the higher energy electron recoil calibrations: 



51 — [Snr / See) ^es{Enr) Ly Enr 



(1) 
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FIG. 1: £cff as a function of recoil energy. The points correspond to the measurements of Manzur 
et al. [42] with statistical and systematic errors in C^s ^ indicated (uncertainties in recoil energy 
not shown). Solid curves show the fiducial C^s dependence used in this work. Filled regions/dashed 
curves indicate the la variation in the C^s dependence. At 3.9 keVnr (the lowest energy data point), 
£eff = 0-0'''3too25io 026 ~ 0.073 ± 0.037. At lower energies three cases are examined: constant £eS 
(blue) at the above value, CeS falling linearly to zero at zero energy (red), and £efr equal to zero 
(green). Above 3.9 keVnr, the gray curve and region are used in all cases. Note linear relationships 
appear curved in the figure due to the logarithmic scaling. 



Here, Ly is the light yield in PE/keVee for 122 keVee 7-rays^. CesiEnr) is the scintillation 
efficiency of nuclear recoils relative to 122 keVee 7-rays in zero electric field; this factor is a 
function of the nuclear recoil energy. Since there is an applied electric field in the experiment, 
which reduces the scintillation yield by quickly removing charged particles from the original 
interaction region, two additional factors must be taken into account: See and Snr are the 
suppression in the scintillation yield for electronic and nuclear recoils, respectively, due to 
the presence of the electric field in the detector volume. The quantities See, Snr, and Ly are 
detector dependent; Ces is not. 

Recent comments have drawn attention to the role that Ces determinations play in setting 
experimental constraints for Xenon-based detectors (32l-[3^ . A variety of CeS measurements 
have been made over the years but limited statistics and systematics issues have so 

far prevented a clear picture from emerging as to the behavior of C^s at low recoil energies. 
There are two primary issues in debate: (1) Which of the C^s measurements should be used 
as a basis for analyzing direct detection results? and (2) Measurements of Ces have only been 
made at energies above some minimum; what is the behavior of Ces at low energies, where no 
measurements have as yet been made? For the first issue, the XENONIOO collaboration has 



chosen to use a global fit to multiple Ces measurements in their analysis, whereas Ref. [32 



suggests that the recent measurements by Manzur et al. [42[ should be used; in both cases, 
£efr measurements are based upon fixed energy neutron scatters. We do not contribute 



^ The unit "keVee" refers to the electron-equivalent energy in keV, the amount of energy in an electron 
recoil event that would produce a given scintillation signal in the detector (whether or not the scintillation 
was, in fact, produced by an electron recoil). The unit "keVnr" refers to the nuclear recoil energy in keV. 
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to the debate as to which Ccs data sets are most appropriate; however, in the interest of 
examining the most conservative XENON constraints, we use the Manzur et al. data alone 
in our analyses. 

The choice of £eff measurements to use in the XENON analyses has a significant impact 
on the resulting constraints for low WIMP masses. The Manzur et al. data yield the lowest 
values for £eff among the fixed-energy neutron scatter measurements, implying the highest 
recoil energy thresholds and therefore the lowest sensitivity for the XENON detectors to low 
mass WIMPs (which generate only low energy recoils). The Manzur et al. data is shown in 
Figure [H 

A comment is in order about the lower ZEPLIN-III £eff measurement 43[ represented 
as a band in Fig. 1 of j32|. ZEPLIN fits a nonlinear £eff model to their broad spectrum 
nuclear recoil calibration data to obtain £eff curves that were used in their analysis. These 
fits suggest a constant £eff at recoil energies above ~30 keVnr, with £efr sharply falling at 
energies below ~20 keVnr and approaching zero at ~7-8 keVnr; see Figure 15 of Ref. ^sf and 
the accompanying text. Thus the suggestion has been made by nonmembers of the ZEPLIN 



team |31| that £eff should be taken to be zero below ~8 keVnr as a conservative model of 
LeS- This £eff model would yield significantly weaker XENON constraints relative to what 
we have referred to as conservative models based on the Manzur et al. measurements of £eff 



42(1 . However, the members of the ZEPLIN experiment themselves do not advocate their fits 



as being an indicator of C^s behavior at recoil energies below ~ 8 keVnr j4J|. In addition, 
the dependence of these curves on statistical and systematic uncertainties has not been 
fully determined, where these uncertainties can significantly impact the lowest recoil energy 
portion of their fits. The ZEPLIN-III dark matter analysis is, in fact, mainly insensitive to 
the low recoil energy portion of their £efr curves. Further discussion of the ZEPLIN data 
can be found in the Appendix. Moreover, as explained in detail in the Appendix, below the 
recoil energies of 7 keVnr, the Manzur et al. measurements are incompatible with £efr = 
at far more than the Sa level. We consider the Manzur et al. measurements more reliable 
than the ZEPLIN-III estimate of £efr at low energies. 

The second issue in debate is how £efr behaves at energies below where measurements have 
been made. Most £efT measurements are at recoil energies above 5 keVnr; Manzur et al. have 
a measurement at 3.9 ± 0.9 keVnr. The £efr behavior below these energies is unclear from 
an experimental and theoretical standpoint, at least at the precision necessary for use in a 
WIMP constraint analysis. The XENON collaboration has suggested that C^s measurements 
are consistent with £efr being effectively constant at low recoil energies, at least at energies 



where recoils may contribute to their signal |12| . |40| . |4l| . Various £eff measurements are also 



consistent with an Ces that decreases as one goes to lower recoil energies; see e.g. Sect. V of 
Ref. 42j which provides a theoretically motivated empirical model of such a decreasing £eff- 
Furthermore, Ref. [32| states that "the mechanisms behind the generation of any significant 
amount of scintillation are still unknown and may simply be absent at the few keVnr level." 
Given this uncertainty we use three different extrapolations of £efr at low energies: constant, 
decreasing as one goes to lower recoil energies, or just zero. 

We choose as our fiducial £efr model a piecewise linear interpolation between the central 
Manzur et al. values at their measured energies, shown in Figure [H In addition, we will also 
examine similarly constructed £eff models using the la uncertainties in the Manzur et al. 
measurements^. The choice of linear interpolation vs. a quadratic interpolation or spline 



The statistical and systematic errors in the measured LcS are added in quadrature with the upper and 
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FIG. 2: The average 51 signal as a function of nuclear recoil energy for XENONIO (left) and 
XENONIOO (right). Curves and regions correspond to the C^s models shown in Figure [TJ 

fit to the £eff points has a negligle impact on the generated constraints compared to that 
from the la variations in the Ces points themselves. Below recoil energies of 3.9 keVnr, 
the lowest Manzur et al. measurement^, we examine three behaviors for £efT, also shown in 
Figured] (1) a constant ^efr, (2) an C^s that goes linearly to zero at zero recoil energy, and 
(3) an £eff that is strictly zero. Even if the scintillation goes to zero at some low but finite 
recoil energy, there is no reason to expect this to occur above ~2-3 keVnr; the measurements 
of £eff provide no indication of an abrupt (rather than gradual) falling of CeS at energies 
just below where the measurements exist. As such, the third case is perhaps unrealistically 
conservative, but never-the-less provides the most conservative case. In addition, the use 
of this case will allow us to examine the contribution of low energy recoils in generating 
constraints. The average 51 signals as a function of the nuclear recoil energy Enr that 
correspond to these Ces models are shown in Figure [2] for XENONIO and XENONIOO. 

In the interest of examining the most conservative XENON constraints, we base all three 
cases on the data from Manzur et al. [i^ . Of the existing data sets, the Manzur et al. data 
yield the lowest values for Ces, implying higher recoil energy thresholds for the XENON 
experiments, and thereby reducing the sensitivity of XENON to low mass WIMPs. 

B. Channeling Effects in DAM A/LIBRA 

The channeling effect is of crucial importance when considering the compatibility of 
DAMA with other experimental results as this effect has the potential to significantly alter 
the WIMP masses and cross-sections which are compatible with the DAMA modulation 
signal. 



lower uncertainties averaged; uncertainties in the corresponding recoil energies for those measurements 
have been neglected. 

Manzur et al. measure £cff = 0.073j:!^;J^^5lJ];226 ~ 0-073 ± 0.037 at 3.9 keVnr, where the two errors in the 
first case are the statistical and systematic uncertainties, respectively, and the second case is the combined 
uncertainty as described in the previous footnote. 
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FIG. 3: Upper bounds to the channeling fraction at a temperature of 293 K for Na (sohd Hues) 
and I (dashed lines) recoiling ions in a Nal crystal for two different models of the temperature 
effect in the lattice parameterized with c = 1 (black) and c = 2 (green or gray). No dechanneling 
processes are taken here into account. To be conservative, in this paper we will use the c = 1 results 
presented here, as they yield the largest change in the DAMA compatible re gion s of parameter 
space relative to the no-channeling case. This figure is reproduced from Ref. 13(], where further 
details may be found. 



Channeling and blocking effects in crystals refer to the orientation dependence of charged 
ion penetration in crystals. In the "channeling effect," ions incident upon a crystal along 
symmetry axes and planes suffer a series of small-angle scattering that maintain them in 
the open "channels" in between the rows or planes of lattice atoms and thus penetrate much 
further into the crystal than in other directions. Channeled incident ions do not get close 
to lattice sites, where they would be deflected at large angles, and they lose energy almost 
exclusively into electrons. The "blocking effect" consists in a reduction of the flux of ions 
originating in lattice sites along symmetry axes and planes, creating what is called a "block- 
ing dip" in the flux of ions exiting from a thin enough crystal as a function of the exit angle 
with respect to a particular symmetry axis or plane. The potential importance of the chan- 
neling effect for direct dark matter detection was first pointed out by H. Sekiya et al. fZsj 
and subsequently for Nal(Tl) by Drobyshevski ^46i] and by the DAMA collaboration |47i] . 
When Na or I ions recoiling after a collision with a dark matter WIMP are channeled, their 
quenching factor^ is approximately Q = 1 instead of Q/ = 0.09 and QNa = 0.3, since they 
give their energy to electrons. The DAMA collaboration |47[ estimated the fraction of chan- 
neled recoils and found it to be large for low recoiling energies in the keV range. Using this 
evaluation of the channeling fraction, the regions in cross-section versus mass of acceptable 



^ The quenching factor Q is the ratio of ionization or scintillation produced by a nuclear recoil event in a 
crystal relative to that produced in an electron recoil event of the same energy. This is analogous to the 
CcS factor in liquid Xenon and is likewise used to reconstruct the nuclear recoil energy from the observed 
ionization/scintillation of an event. 
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WIMP models in agreement with the DAMA data were found to be considerably shifted 
towards lower WIMP masses and cross-sections. 

However, the DAMA calculation of the channeling fraction did not take into account 
that the recoiling lattice ions start initially from lattice sites (or very close to thernWind, 
therefore, blocking effects are important. In fact, as argued originally by Lindhard [4a], in 
a perfect lattice and in the absence of energy-loss processes, the probability of a particle 
starting from a lattice site to be channeled would be zero. The argument uses statistical 
mechanics in which the probability of particle paths related by time-reversal is the same. In 
a perfectly rigid lattice, the fraction of channeled recoils would, in fact, be zero. However, 
the atoms in a crystal are actually vibrating about their equilibrium positions in the lattice. 
It is this displacement from equilibrium that allows for a non-zero channeling probability of 
recoiling ions. The vibration amplitude increases with the temperature, thus the effect is 
temperature dependent: in general the channeling fraction increases with temperature. 

Upper bounds to the recoiling channeling fractions in Nal(Tl) crystals at 20°C were 
obtained in Ref. [isj, using analytic models of channeling developed since the 1960's, when 
channeling was discovered (see for example Refs. 48N50| and references therein). These 
upper bounds on the channeling fractions were obtained with temperature effects taken into 
account not only through the vibrations of the colliding nucleus but also in the lattice. The 
latter depend on the parameter c (see Ref. 13| for details) which in the relevant literature 
is found to be a number between 1 and 2, with 1 giving the largest channeling fractions (see 
Figure [3l reproduced from Ref. fisf). 

The fractions shown in Figure E] are also an upper bound in that no dechanneling mecha- 
nism has been taken into account to compute them. The collisions with Tl impurities would 



take channeled ions out of their channel, and this process is not included (see Ref. [13| for 
further explanations). 



III. RESULTS AND DISCUSSION 

In Figure m we show the WIMP masses and SI cross-sections compatible with the DAMA 
modulation signal both with and without channeling included; contours are shown for regions 
compatible at the Tcr, 5a, 3a, and 90% level (in order from larger to smaller regions). For 
the channeling, we use the largest channeling fractions shown in Figure [3] as they provide 
the largest potential effect on the DAMA constraints. Figure H] shows that even in this case 
there is negligible difference between the channeling and non-channeling scenarios except 
for regions incompatible with DAMA at greater than the 5a level. Even in these cases, the 
difference lies only at WIMP masses below 4 GeV and at relatively high SI cross-sections. 
As channeling is a negligible effect, we do not further include it. 

Compared to our previous analysis in 0] , the current study takes advantage of additional 
recently released DAMA data. The effect of the additional data has been to sharpen the 
regions in parameter space that match the data. For example, at 5a, there are now two 
completely separate regions (peaked at different WIMP masses) that were previously joined. 
In our current work we also display a 7a contour in which the two regions are again connected. 
We remind the reader that we are using the goodness-of-fit statistic described in detail in 
Ref. 0. 

Our main results are shown in Figs. corresponding to the three cases for the 

behavior of £efr at low recoil energies. The solid gray contours indicate the WIMP parameters 
compatible with the DAMA modulation within the 5a, 3a, and 90% level; the 5a DAMA 
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FIG. 4: WIMP masses and spin-independent (SI) cross-sections compatible with the DAMA mod- 
ulation signal and total number of events, determined with (dashed green) and without (solid 
orange) the channeling effect included. The largest channeling fractions shown in Figure [3] (taken 
from Ref. 13|]) are used here for the channeling case. Comparing the cases with or without chan- 
neling, we find negligible difference in the DAMA modulation regions at the 90%, 3a, and 5a 
levels; only the 7a contours differ and only for WIMP masses below 4 GeV. The lower and higher 
mass DAMA regions correspond to parameters where the modulation signals arise from scattering 
predominantly off of Na and I, respectively. 



region is also shaded light gray. The (filled) pink contour corresponds to the 7-12 GeV WIMP 
mass region suggested by CoGeNT (we reiterate that we have not reanalyzed their data and 
simply display their published region here). CDMS, DAMA (total events), XENONIO, 
and XENONIOO curves indicate regions for which the WIMP parameters are excluded at 
the 90% level (the parameters above these curves are excluded). The solid green region for 
XENONIO and solid purple region for XENONIOO do not indicate regions compatible within 
a given level (as opposed to the DAMA and CoGeNT regions); they instead indicate how 
the 90% exclusion constraints vary with the la level uncertainties in the Ces measurements. 
Overlapping XENONIO and XENONIOO la regions are shown in blue. 

For the fiducial (central value) Ces model in the case where it is constant below 3.9 keVnr, 
shown in Figure O the XENONIOO constraint excludes all of the DAMA 3a region, but only 
the portion of the CoGeNT region with WIMP masses above 9 GeV. We note that, because 



we use only the Manzur et al. C^s data |42[, this constraint is weaker than that presented 
by XENONIOO [l2|. If the la uncertainties in C^s are included, XENONIOO could exclude 
nearly all of the DAMA 5cr region and the entire CoGeNT region, if the largest value of £efr in 
the la region is taken. On the other hand, it might exclude only the CoGeNT region above 
11 GeV and not even all of the DAMA 90% region, if the lowest value of £eff in the la region 
is used. However, the CDMS constraint, unaffected by the issues with £eff, constrains the 
same CoGeNT region as the fiducial XENONIOO case here, with a slighter weaker constraint 
on the DAMA region (incompatible with the DAMA 2a region, not shown). 
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FIG. 5: XENONIO (green) and XENONIOO (purple) 90% C.L. constraints for a constant £efr at 
recoil energies below 3.9 keVnr. The solid curves are the constraints using the central values of >CefF 
as described in the text; dashed curves and lighter filled regions indicate how these 90% constraints 
vary with the la uncertainties in C^s- The blue region indicates an overlap between the XENONIO 
(green) and XENONIOO (purple) la regions. Also shown are the CDMS constraint (orange curve), 
DAMA modulation compatible regions (gray contours/region), and the CoGeNT 7-12 GeV region 
(pink contour/region). The lower and higher mass DAMA regions correspond to parameters where 
the modulation signals arise from scattering predominantly off of Na and I, respectively. 
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FIG. 6: Same as Figure[5l but taking C^s to fall linearly to zero for recoil energies below 3.9 keVnr. 
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FIG. 7: Same as Figure O but taking C^s to be zero for recoil energies below 3.9 keVnr. 

For the fiducial (central value) C^s model in the case where it falls linearly to zero at zero 
recoil energy, shown in Figure |6l the XENONIOO constraint again excludes nearly all of the 
DAMA 3(7 region and the portion of the CoGeNT region with WIMP masses above 9 GeV. 
The la variations in the C^s measurements also yield a similar variation in the XENONIOO 
constraint as they did in the constant Ces case. The most extreme case, taking Ces to be 
zero below 3.9 keVnr, yields similar XENONIOO constraints as the other two cases, as seen 
in Figure [TJ although the constraint using the la upper values of CeS does not quite exclude 
the full CoGeNT region, leaving a narrow window at WIMP masses of 7-8 GeV. It should 
be emphasized, however, that the linearly falling Ces case is already conservative and taking 
CeS to be zero below 3.9 keVnr is perhaps unrealistically conservative. 

The XENONIOO constraints are nearly identical in the DAMA and CoGeNT regions for 
all three cases of low energy Ces behavior. In fact, the constraints based on the central 
and la lower values of CeS are identical; only when using the upper la CeS values do the 
constraints differ. There are two main reasons for the similarity among the constraints: (1) 
the imposed (SI) > 1.0 PE cutoff and (2) the small potential contribution from recoil events 
with energies below 3.9 keVnr where the C^s models differ. As can be seen in Figure |21 
a recoil energy of 3.9 keVnr yields an average 5*1 signal of 1.0 PE in XENONIOO when 
using any of the three fiducial Ccs models. With the (5*1) cutoff, there is no contribution 
from recoils at energies below 3.9 keVnr where the fiducial Ces models differ; thus, these 
constraints are identical. When using the la lower values of Ces, no recoils below 5.9 keVnr 
are included, so the lesser constraining portion of the la XENONIOO constraint bands shown 
in the figures are likewise identical. On the other hand, when using the la upper values 
of CeS, the (5*1) > 1.0 PE cutoff corresponds to recoil energies of 2.5, 3.1, and 3.9 keVnr 
for the constant, linearly falling, and zero low energy Ces models, respectively. In this case, 
low energy recoils contribute to the constraints. However, these low energy recoils can make 
only a small contribution to the observed signal, as will be discussed below. 

As the potential effect of the low energy Ces behavior on the XENONIOO constraints is 
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FIG. 8: XENONIO and XENONIOO constraints when relaxing the {SI) > 1 PE cutoff. Only the 
constraints for the fiducial (central) £efr niodels are shown; from left to right, these correspond to 
constant, linearly falling, and zero CeS below 3.9 keVnr. To illustrate the potential effect of the low 
recoil energy behavior of C^s, these contraints are generated by arbitrarily assuming the nuclear 
recoil band cut efficiency is constant at low recoil energies. This efficiency should actually fall 
at sufficiently low recoil energies (for XENONIO, at energies that yield an average SI somewhere 
below 1 PE); these constraints should therefore not be taken as true constraints on the WIMP 
mass and cross-section. The actual constraints based upon a proper accounting of the efficiencies 
at low recoil energies will fall somewhere between the constraints shown here and in the previous 
figures. 

masked by the (5*1) > 1.0 PE cutoff, we show in Figure [H] the XENONIOO constraints for 
the three fiducial Ces niodels when this cutoff is relaxed. In this figure, we have arbitrarily 
assumed the nuclear recoil band cut efficiency is constant at low recoil energies. In reality, 
this efficiency should fall at very low recoil energies and this approximation becomes inap- 
propriate at recoil energies that yield (5*1) somewhere below 1 PE. For this reason, these 
constraints should not be taken to be valid constraints; we show them only to illustrate 
the potential effect of low energy recoils and the low energy Ces behavior. With all cut 
efficiencies properly taken into account, the true constraints would lie somewhere between 
the constraints shown in this figure and those shown in the previous figures. 

With the relaxing of the (SI) cutoff, the XENONIOO constraints are still nearly identical 
in the DAMA and CoGeNT regions for all three cases of low energy Ces behavior. These 
constraints are very similar to the ones found in the previous figures, when a (SI) cutoff was 
included, and are actually identical for the zero Ces model as this model has no contributions 
from recoils with (SI) < 1.0 PE anyways. The three cases only begin to differ significantly 
in the low mass, high cross-section parameter space located around and above the DAMA 
regions in this figure. This can be explained by the XENONIOO SI analysis threshold of 
4 PE's (the full analysis range is 4-20 PE's). In the absence of a finite energy resolution, 
this corresponds to a nuclear recoil energy of 9.5 keVnr in our fiducial C^s models, well 
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into the energy range where the Ces behavior is known. With a Poisson fluctuation in 
the number of observed PE's, recoils at lower energies have a finite chance of producing 4 
or more PE's and falling into the analysis range, even if the average number of PE's for 
events at those energies is below 4. However, at 3.9 keVnr, the average expected number 
of PE's is 1.0; only 1.9% of such events yield 4 or more PE's. Recoils of 3 keVnr yield an 
average number of expected PE's of 0.79 and 0.61 for the constant and falling Ces cases, 
respectively, with corresponding probabilities of being observed (4+ PE's) of 0.87% and 
0.36% (the third case, zero Ces, produces no PE's at these energies). The small fraction 
of recoil events with energies below 3.9 keVnr that will be observed in the analysis range 
means that their contribution is only significant when there are essentially no events at 
higher energies (due to low WIMP masses and a finite escape velocity in the halo), but to 
produce a sufficient number of events to fall into the analysis range requires a very large 
number of WIMP scatters in the ~l-4 keVnr range, which requires a high WIMP cross- 
section. Thus, even when using an overly optimistic nuclear recoil band efficiency, the three 
£eff cases can only result in different constraints in the low mass, high cross-section region. 
This is not necessarily the case for C^s curves based on measurements that yield values 
higher than Manzur et al, as this would push the analysis range corresponding to 4-20 PE's 
to lower recoil energies; such cases, however, inevitably move the XENONIOO constraints 
to the left. In any case, the most significant issue in the XENONIOO analysis is the choice 
of Ces measurements used to determine the Ces dependence, not the Ces behavior at low 
energies. 

We now turn to the XENONIO bounds. We have reanalyzed the XENONIO results in 
terms of the same CeS models as used for XENONIOO and discussed in the previous section; 
our results differ from those shown by the XENON collaboration due to the difference in 
Ces used in their analyses and ours. The XENONIO results are important because of the 
lower 5*1 threshold of about 2 PE's used in that analysis, which corresponds to 4.6 keVnr 
nuclear recoil energies in our fiducial Ces models (neglecting Poisson fluctuations), much 
lower than the 9.5 keVnr of the XENONIOO 4 PE threshold. Because of the lower threshold, 
the behavior of Ccs at low recoil energies is relevant in producing the XENONIO constraints 
as Poisson fluctuations allow for a non-trivial probability of seeing 2+ PE's for recoil energies 
below 3.9 keVnr. 

For the constant £efr case shown in Figure El the lower threshold allows for a stronger 
sensitivity to lower WIMP masses for XENONIO relative to XENONIOO. The fiducial case 
excludes at the 90% CL. all of the CoGeNT region and DAMA to the 5a contour. When 
the la uncertainties in the Ces measurements are taken into account, the constraints relax: 
DAMA is excluded to only the 3a contour and the some of CoGeNT region at WIMP 
masses below 9 GeV survive. The XENONIO constraints mildly weaken if Ces is taken to 
fall linearly to zero below 3.9 keVnr, as seen in Figure [61 The fiducial case still excludes 
all of the CoGeNT region and DAMA to about the Aa contour (not shown). The Ces lo" 
band here allows the same DAMA and CoGeNT regions to survive as with the constant Ces 
case. For the case where C^s is zero below 3.9 keVnr, shown in Figure [71 the XENONIO 
constraints further weaken and approach the XENONIOO constraints as the low energy 
events are essentially turned off and the lower XENONIO threshold becomes less relevant. 
Again, we note that this last case (zero CeS at low recoil energies) is an extremely conservative 
case. 

As with XENONIOO, the potential effect of the low energy Ces behavior on the XENONIO 
constraints is limited by the imposed (5*1) > 1.0 PE cutoff. We also show in Figure [H the 
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XENONIO constraints for the three fiducial £eff niodels when this cutoff is relaxed. The 
same caveats apply: the nuclear recoil band cut efficiency that is used is not appropriate for 
the full recoil energy range that it is applied over, so these do not represent valid constraints. 
Again, these constraints are only used to illustrate the potential impact of the low energy 
£eff behavior on XENONIO constraints. The actual constraints when all efficiencies are 
accounted for properly would fall somewhere between the constraints shown in Figure [H] and 
those shown in Figs. ©-(ITj). 

With the (5*1) cutoff relaxed. Figure [8] shows how the low threshold allows for a strong 
XENONIO sensitivity to lower WIMP masses. This is particularly evident with the constant 
£cff case where, as can be seen in Figure [2l recoils of energy 1 keVnr yield an average SI 
signal of 0.4 PE; ~6% of such recoils will produce the necessary 2+ PE. For the falling Ccs 
case, that same average 5*1 signal of 0.4 PE occurs at a higher recoil energy of 2 keVnr, but 
this energy is still sufficiently low to provide sensitivity to low mass WIMPs. The presence of 
these non-trivial Poisson ffuctuations at low recoil energies leads to a very strong dependence 
of the XENONIO constraints on the low energy Cgg behavior. This should remain the case 
even when the various efficiencies are handled properly, though not quite to the degree 
shown in Figure [HI In particular, when the proper efficiencies are included, the XENONIO 
constraints in the constant and falling Ces cases should gain an upward curve at low WIMP 
masses, as seen with the other constraints, rather than the current linear appearance. These 
linear portions of the constraints at low WIMP masses (as they appear with the logarithmic 
scaling of the figure) continue to arbitrarily low WIMP masses; however, they arise from the 
Poisson tails of increasingly smaller energy events that would be suppressed when using the 
proper efficiencies. 

Though we have not included it in this work, CDMS Silicon data may provide further 
constraints on the DAMA and CoGeNT regions and should be considered in a full discussion 
of compatibility between the various experimental results. The CDMS Silicon results will 
be considered in future work. 

In summary, we have examined a number of subtleties relevant to direct detection studies 
of low mass WIMPs. In the interest of examining the most conservative XENON constraints, 
we have used the Manzur et al. {i^ data alone in our analyses. Of the existing data sets, the 
Manzur et al. data yield the lowest values for £efr; implying higher recoil energy thresholds 
for the XENON experiments, and thereby reducing the sensitivity of XENON to low mass 
WIMPs. We find that, when basing the £eff curves on these Manzur et al. measurements, the 
behavior of £eff at low energies (less than 3.9 keVnr) has negligible effect on the XENONIOO 
constraints in the regions of interest for DAMA and/or CoGeNT. For XENONIOO, the choice 
of data sets upon which the £efr dependence is based is more important than the extrapolated 
behavior of £cfr at low recoil energies. The strongest bounds are from XENONIO, rather 
than XENONIOO, due to the lower energy threshold. For reasonable choices of C^s and for 
the case of spin independent elastic scattering, we find that XENONIO is incompatible with 
the DAMA/LIBRA Sa region and severely constrains the CoGeNT 7-12 GeV WIMP mass 
region. 
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Appendix: Comment about our choice of relevant parameters and response to 

critique 

Shortly after the first version of our paper was released, Ref. [31] appeared commenting 
on it. We include here a detailed explanation of the XENONIO efficiencies and cuts we are 
using, as well of our use of the Manzur et al. measurements as a conservative choice for 
£eff (instead of the ZEPLIN-III measurements) which we believe are relevant in view of the 
comments expressed in Ref. (3l| . 



1. XENONIO efficiencies and cuts 

The 51 peak finding efficiency factor rjsi had not been included in the XENONIO analysis 



in the first version of our paper (as correctly pointed out in Ref. 3l|) and has been accounted 
for in this revised version. However, this leads to only a moderate weakening of the XENON 
constraints by shifting the bound on the cross section upward by less than a factor of two. It 
is important to point out that including this effect does not weaken the constraints upwards 
in cross-section by 2-3 orders of magnitude, contrary to the claim in the critique given in 
Ref. |3l|. One can understand the small magnitude of the effect with the following reasoning. 

A valid signal event in the XENON detectors is required to produce coincident scintillation 
in at least two PMTs. The rjsi factor accounts for experimental limitations in identifying and 
reconstructing at least two PMT contributions to the overall 5*1 signal of a recoil event. The 
SI signal is determined from the area under the peaks produced in the electronic readout 
of the PMTs. Due to digitization of the signal and intrinsic PMT performance, the size 
and shape of the peaks will vary^; see Fig. 14 of Ref. (sij for an example. Small or poorly 
shaped peaks may fail to be properly tagged as a PE peak in a PMT. Using the more 
conservative estimate found in Ref. j25|, only rfsi ~ 60% of 2 PE events will have both PE 
peaks properly tagged and identified as coincident. As the number of PE's increases, the 
probability of passing the two-fold PMT requirement rapidly rises to 100%. 

The significance of the rjsi factor on XENONIO constraints can be easily estimated. If 
this factor is conservatively assumed to be 60% over the entire 5*1 analysis range (2+ PE) 
instead of just at low 5*1, the expected number of events passing the XENONIO cuts at all 
energies should fall to 60% of the number of events expected without rjsi applied. As the 



^ The measured 51 signal from a single PE as determined from the area of these peaks is 1.0 ± 0.6 PE. 
Thus, SI can take on non-integer values even though it is given in terms of number of PE. 
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FIG. 9: Impact of (i) the SI peak finding efficiency rjsi (dashed Unes) and (ii) the {SI) > 1 cutoff 
(dotted hnes) on the XENONIO constraints. Only the constraints for the fiducial (central) C^s 
models are shown; from left to right, these correspond to constant, linearly falling, and zero C^s 
below 3.9 keVnr. One can see that the second effect is by far the stronger one. 



recoil rate at a particular WIMP mass is proportional to the scattering cross-section, the 
original recoil spectrum can be exactly reproduced by shifting the cross-section upwards by 
a factor of 1/0.6 ~ 1.7. Thus, if a WIMP mass and cross-section is excluded by XENONIO 
at some given CL without the rjsi factor applied, a WIMP with the same mass and xl.7 
higher cross-section would yield the same excluded spectrum with the rjsi factor included. 
In this case, the XENONIO exclusion curves shift upwards in cross-section by xl.7. As t]si 
is higher than 60% for 5*1 larger than 2 PE, the shift in cross-section will be even smaller 
at heavier WIMPs where a significant number of high 5*1 events are expected. 

We demonstrate in Figure |9] the weakening of the XENONIO constraints when rjsi is 
included. We show only the three fiducial (central) C^s niodels described in our paper; 
constraints without and with rjsi included are given by solid and dashed curves, respectively. 
For a given WIMP mass, the constraint increases by ~ xl.7 in the scattering cross-section 
at low WIMP masses and by less at higher WIMP masses. This is a relatively modest 
weakening of the constraint given the logarithmic scaling of the exclusion curve figures. 



Please recall that we have used the more conservative estimate of rjsi found in Ref. |25 
Had we used the risi factor given in Ref. [4o[, the change would have been even milder. 

We will clarify an issue that is perhaps the source of what appears to be an erroneous 
application of the rjsi factor in Ref. [3l| that yielded a greater weakening in the XENONIO 
results than expected. Given a recoil at some energy Enr, Eqn. ([T]) gives the average expected 
5*1 scintillation signal (5*1). Given the measured 5*1 in an event, Eqn. ([T]) can be inverted 
to obtain a reconstructed recoil energy E'. Due to the discreteness of the produced PE's 
and the variation in the 5*1 peaks,the observed prompt scintillation signal 5*1 is a random 
value whose expectation value is (5'1). Thus, the reconstructed recoil energy E' is also 
a random value and provides only an estimate of the (unknown) true recoil energy, Enr- 
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The distinction between E' and Enr is important. With the PE fluctuations and 5*1 peak 
variations, the true recoil energy cannot be precisely determined for any given recoil event. 
E' is simply an estimate of the likely recoil energy that produced an event and is understood 
to have some intrinsic uncertainty attached to it. At high recoil energies, E' Enr and £" 
can be taken as a good approximation of the actual recoil energy. At low recoil energies, E' 
may differ signiflcantly from Enr on an event by event basis. 

The rjsi factor is sometimes given as a function of recoil energy. However, as rjsi is really 
a function of the 5*1 signal, this recoil energy refers to the reconstructed energy E' of an 
event (recall E' and 5*1 have a 1:1 mapping through Eqn. ([1])) and does not refer to the 
actual recoil energy Enr- Neglecting the 5*1 peak variations and efficiencies other than rjsi, 
the fraction of events / that exceed the 2 PE threshold from scatters at a recoil energy that 
yields an average 5*1 signal (5*1) is given by: 

f{{Sl)) = vsi{2)P{2\{Sl)) + r/5i(3)P(3|(51)) + vsi{^)P{A\{Sl)) . . . , (A.l) 

where rjsi is taken as a function of 5*1 and P(A;|(5'1)) is the Poisson probability of seeing 
k PE with an average of (SI). If the risi factor were assumed to be a function of the actual 
recoil energy Enr, then 

f{{Sl)) = VsiiiSl)) [P{2\{S1)) + P{3\{S1)) + P{4\{S1)) . . .] (A.2) 

would erroneously be taken as the fraction of events exceeding the 5*1 threshold. Use of 
Eqn. ( 1A.2|) in place of Eqn. ( ]A.1|) will improperly yield greatly weakened XENONIO con- 



straints at low WIMP masses. 

A second change has been made in the present version of this paper. In the first version 



we assumed the same nuclear recoil band efficiency given for 2 < 5*1 < 5 in Ref. [11 
applied to all Poisson fluctuated events that appeared in that range. This was an optimistic 
assumption that may overestimate the XENONIO sensitivity at low WIMP masses. There 
are two data cuts relevant for low energy recoils that yield 2 < 5*1 < 5: a recoil event must 
have 52 > 300 PE (52 threshold) and 1.88 < \og^o{S2/Sl) < 2.40 (nuclear recoil band) to 
be accepted as a valid event. The latter requirement is a cut designed to exclude electron 
recoil background events which tend to produce higher values of 5*2/5*1 than nuclear recoils; 
the range of values here accept the lower ~45-50% of the distribution of 52/51 expected for 
nuclear recoils, as determined from calibration data. At very low recoil energies, an event 
that produces an upward fluctuated 51 > 2 PE and would otherwise be considered a valid 
event might fail to produce enough 52 signal to pass the 52 threshold. The efficiency can 
thus be expected to fall at low recoil energies. To avoid problems with these two data cuts, 
we have added a cutoff to our analysis and ignore recoil energies for which (51) < 1.0 PE. 
We also take the same 47% nuclear recoil band efficiency in the 2-5 PE 51 bin for all events 
in this 51 range stemming from low energy recoils. The addition of this cutoff significantly 
weakens the XENONIO constraints for WIMP masses below ~6 GeV for the constant and 
falling Ces models discussed in the paper, as seen by the dotted curves in Figure [E This 
cutoff has a far greater impact on the constraints than the inclusion of the risi factor and 
is the dominant source of the change in XENONIO constraints from the first version of this 
paper. The third £eff model, zero below recoil energies of 3.9 keVnr, is unaffected by this 
cutoff as it does not yield any observable events at low energy anyways. 

The (51) > 1.0 PE cutoff and 47% assumed efficiency we have adopted here are con- 
servative for two reasons: (1) The 51 and 52 fluctuations are independent. This means 
that the Poisson fluctuated events that yield higher 51 than the average will have lower 
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S2/S1 ratios on average and are more likely to pass the nuclear recoil band cut. The up- 
ward Poisson fluctuated events for 1.0 < (5*1) < 2.0 can survive this cut as much as 70-80% 
of the time, higher than the assumed 47%. (2) Our choice of a cutoff at (SI) = 1 PE is 
due to our limited ability to examine the efficiencies at lower recoil energies, not due to an 
expected lack of events at these low energies. We note that there will also be fluctuations 
in the number of ionization electrons that lead to the 5*2 signal even at very low recoil 
energies, there may be a non-zero probability of both 5*1 and 5*2 exceeding their respective 
thresholds and producing events in the analysis region that pass all cuts. For the constant 
CeS case at a WIMP mass of 5 GeV and scattering cross-section of 10~^ pb, a point not 
excluded by the conservative bounds in Figure [9l about 20,000 recoils would be expected for 
recoil energies corresponding to 0.5 < (5*1) < 1.0. Even with a very small efficiency for such 
recoils, some number of these events would be expected to have fluctuations that put them 
in the XENONIO analysis region. If the nuclear recoil band efficiency were only ~1% over 
this range, 10-|- events would be expected to pass the various cuts; the lack of such events 
in the data would rule out these WIMP parameters. 



2. ZEPLIN-III C^ff models 



Refs. [32| and [3l[ suggest using the the ZEPLIN-III C^s measurement [43|, represented 



as a band in Fig. 1 of Ref . [32| , as a conservative choice of Ces- Thus, we discuss here this 
possibility. 

ZEPLIN fits a nonlinear Ces niodel to their broad spectrum nuclear recoil calibration 
data to obtain CeS curves that were used in their analysis. These flts suggest a constant CeS 
at recoil energies above ~30 keVnr, with Ces sharply falling at energies below ~20 keVnr 
and approaching zero at ~7-8 keVnr; see Figure 15 of Ref. |43j and the accompanying text. 



The suggestion has been made [31[ that Ces should be taken to be zero below ~8 keVnr 
as a conservative model of Ces- This Ces model would yield signiflcantly weaker XENON 
constraints relative to what we have referred to as conservative models based on the Manzur 
et al. measurements of £eff (42| . 

We have several reservations about using this ZEPLIN inspired £eff niodel. First, ZEPLIN 
does not provide estimates of Ces below ~7 keVnr and, in fact, states they are limited in 
constraining Ces at these low energies (in the caption of their Fig 15, it is said that "The 
constraints become very weak outside the energy ranges shown."). Second, neither the 
technical details of ZEPLIN's curve-fltting nor estimates of the statistical and systematic 
uncertainties are provided in their paper. Without these, it is unclear with which degree 
of certainty the low energy (~7-10 keVnr) end of the ZEPLIN-III curves should be treated. 
Several sources of error may contribute to the level of uncertainty in these C^s estimates. One 
is the errors associated with the Monte Carlo used to compare with data, which are difficult 
to quantify since they can arise from inaccuracies in the inelastic neutron scattering data 
used in the Geant4 code and obtained from an international library of such data. Another 
issue are systematic errors associated with uncertainties in the position of the neutron source 
near the detector during calibration runs. If the uncertainties are large, then ZEPLIN does 
not have the statistical power to determine the Ces behavior at low recoil energies and 
one should base the C^s curves on measurements that have a better statistical power to 



^ The S2 threshold of ^-^300 PE corresponds to about 12 ionization electrons, each of which produce ^-^25 PE. 
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analyze this Ccs behavior, such as the Manzur et al. data. If the uncertainties are small so 
that the ZEPLIN curve is expected to be an accurate representation of the Ces behavior 
at low energies, then there is a very serious discrepancy between the ZEPLIN and Manzur 
et al. results (see below). The ZEPLIN collaboration claims no leverage on C^s below about 
8 keVnr and do not suggest that it goes to zero in that energy range* [44| • The recoil energy 
interval in the ZEPLIN analysis of 10.7-30.2 keVnr does not include the low energy, low Ccs 
end of their Ces curves. In this analysis region, the ZEPLIN Ces curves are compatible with 
the Manzur et al. measurements within the l-2a level. 



Furthermore, the ZEPLIN inspired model proposed by Ref. [31[, taking C^s = below 
recoil energies of ~8 keVnr, is strongly incompatible with the Manzur et al. measurements 
as well as other fixed energy neutron scattering measurements such as those of Aprile et al. 



41[. We reiterate that our choice to use the Manzur et al. data in our analysis was due 
to the lower C^s values that give more conservative XENON constraints and should not 
be construed as an indication that we regard this data as the most accurate available. We 
make no contribution to the discussion regarding the accuracy of the various low energy 



recoil measurements from neutron beam scattering. See Ref. [52| for a discussion of the 
various potential issues which may affect the Ces determinations in these experiments. In 
any case, much of the discussion below regarding the incompatibility of the Manzur et al. 
data with CeS = applies also to the Aprile et al. measurements. 

Given the Manzur et al. measurement of e.g. CeS = 0.073^o!o25-ao26 (statistical and 
systematic errors, respectively) at a recoil energy of 3.9 ± 0.9 keVnr, one might naively 
conclude this measurement is consistent with C^s = at the ~2-3(T level. However, the 
Ccs value and errors are determined from a fit to the data as a function of £eff- The 
versus Ces curves, such as shown in Fig. 11(c) of Ref. [42}] for a recoil energy of 6 keVnr, are 
not symmetric about the minima (which provide the central values), but rise very rapidly 
as Ces becomes small, effectively diverging at Ces = 0. The result is that the Manzur et al. 
measurements (four of which are below recoil energies of 7 keVnr) are incompatible with Ccs 
= recoil energies of ~8 keVnr, as this would yield values that are imcompatible at far 
higher than the 3a level. We consider the Manzur et al. measurements more reliable than 
the ZEPLIN-III estimate of £eff at low energies. 

To illustrate the incompatibility of the Manzur et al. data with C^s = 0, one can under- 
stand the behavior as follows. Manzur et al. measured Ccs at several recoil energies by 
observing the scintillation response in a Xenon detector for neutrons from a beam of fixed 
energy (2.8 MeV) that scatter in a fixed direction; the recoil energy is fixed by the angle 
of scatter^. If Ces is non-zero at the recoil energy corresponding to a particular scattering 
angle, a histogram of the 5*1 from observed events will generate a peak due to scintillation 
from the single scatter events in the detector. The primary background is double scatter- 
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ZEPLIN used a 9-point spline function to model the curves and a classical maxmium likelihood method 



using a grid scan across a 9 point parameter space to make sure they found a global minimum 44j . 



ZEPLIN only fit their data down to 2 keVee, thus could infer nothing below about 6 keVnr. The middle 



curve in their Fig. 15 43[ gives the "best fit" and the outer two curves indicate regions where "similar 
goodness-of-fit values" were obtained as stated in the caption; the band does not represent a region with 
a particular statistical significance {e.g. 1-a or 90% C.L.) in the uncertainties. 

Due to the geometry of the Xenon detector and width of the scintillator used to detect the scattered 
neutrons, each measurement is sensitive to a finite range of recoil angles and, thus, recoil energies. The 
uncertainties in the recoil energies given for each CcS measurement arise mainly from this range. 
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ing neutrons which can reach the neutron scintillator located at the fixed angle from the 
neutron beam, but scatter with a different energy than expected for single scatter neutrons 
that reach the neutron scintillator; one of the two scatters must occur outside of the active 
volume for the event to be mistaken as a single scatter. These double scattering neutron 
events are expected to occur at an almost negligible rate compared to the single scatter 
signal events. In addition, these double scatter events produce a fairly flat SI distribution. 



not a peak (see e.g. Figure 9 of Ref. [42|]). Figures 11(a) and 11(c) of Ref. 42| show an 
SI histogram and corresponding versus Ces curve, respectively, for recoils centered at 
6 keVnr. The histogram shows a peak of ~700 events (consistent with scintillation from 
single scatter signal events) with a negligible flat contribution (expected from any signiflcant 
double scatter background). The this case is minimized at Ces ~ 0.06. At lower Ces, 

the grows rapidly for this reason: it is difficult to account for this ~700 scintillation 
event peak if the nuclear recoils produce little to no scintillation. Note the systematic issues 
discussed in Manzur et al. generally affect the interpretation of which C^k would produce 
such a peak, but do not provide significant alternate mechanisms for producing such a peak 
aside from the scintillation from low energy recoils (scintillation that can only be produced 
if £eff 7^ 0). Thus, while determining the value of £eff based on the size and shape of such 
5*1 peaks may be limited by statistics or biased by systematic effects, the simple presence 
of such a peak is strong evidence that £eff is non-zero. 

Given that Manzur et al. have observed scintillation peaks at multiple recoil energies 
over ~ 4-8 keVnr, we consider the "conservative" £efr model suggested in Refs. and l3ll . 
with £eff = below ~8 keVnr, to be grossly incompatible with existing data and therefore 
unrealistically conservative. 
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